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We report a new analysis of lattice simulation results for octet baryon masses in 2+1-flavor QCD, 
with an emphasis on a precise determination of the strangeness nucleon sigma term. A controlled 
chiral extrapolation of a recent PACS-CS Collaboration data set yields baryon masses which exhibit 
remarkable agreement both with experimental values at the physical point and with the results of 
independent lattice QCD simulations at unphysical meson masses. Using the Feynman-Hellmann 
relation, we evaluate sigma commutators for all octet baryons. The small statistical uncertainty, 
and considerably smaller model-dependence, allows a significantly more precise determination of 
the pion-nucleon sigma commutator and the strangeness sigma term than hitherto possible, namely 
f-n-jv = 45 ± 6 MeV and a a = 21 ± 6 MeV at the physical point. 
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The light-quark sigma terms provide critical informa- 
tion concerning the nature of explicit chiral symmetry 
breaking in QCD, as well as the decomposition of the 
mass of the nucleon pQ . While these physical observables 
arc difficult to measure with conventional probes, an ac- 
curate knowledge of the sigma terms is of essential impor- 
tance in the interpretation of experimental searches for 
dark matter (2HS]. Dark matter candidates, such as the 
favoured neutralino, a weakly interacting fermion with 
mass of order 100 GeV or more, have interactions with 
hadronic matter which are essentially determined by cou- 
plings to the light and strange quark sigma commutators. 

Experimentally, ct^n is determined from irN scattering 
through a dispersion relation analysis. Traditionally, the 
strange scalar form factor has then been evaluated indi- 
rectly using <j„n and a best-estimate for the non-singlet 
contribution <7o = mi(N\uu + dd — 2ss\N). These tra- 
ditional evaluations have yielded a value for a s as large 
as 300 MeV, compared to 50 MeV for the light quark 
commutator, indicating that as much as one third of the 
nucleon mass might be attributed to non- valence quarks. 
This suggestion appears to be incompatible with widely- 
used constituent quark models, and has generated much 
theoretical interest over the last two decades. 

The traditional method of determination of a s is 
severely limited because it involves the small difference 
between a^N (with its uncertainty) and <7o which is usu- 
ally deduced in terms of SU (3) symmetry breaking. Even 
given a perfect determination of a^N > <?s will have an un- 
certainty of order ~ 90 MeV [7] . For that reason a s has 
been considered notoriously difficult to pin down. In re- 
cent years, the best value for a s has seen an enormous re- 
vision. Advances in lattice QCD have revealed a strange 
sigma term of 20-50 MeV [8HT7] , an order of magnitude 
smaller than was previously believed. 

In this Letter we use the finite-range regularisation 
(FRR) technique to effectively resum the chiral pertur- 
bation theory expansion of the quark mass dependence 



of octet baryons. Fitting the resulting functions to re- 
cent lattice data, we extract the scalar form factors by 
simple differentiation using the Feynman-Hellmann theo- 
rem. Our technique allows comparison with recent direct 
lattice QCD calculations of the flavor-singlet matrix ele- 
ments at unphysical meson masses |19j . with consistent 
results. We report values of a„N — 45 ± 6 MeV and 
a s = 21 ± 6 MeV at the physical point. 

The sigma terms of a baryon B are defined as scalar 
form factors, evaluated in the limit of vanishing momen- 
tum transfer. For each quark flavor q, 

&Bq = niq(B\qq\B); a Bq = a Bq /M B . (1) 

For the nucleon, the so-called ttN sigma commutator and 
the strange sigma commutator are defined by 

ovat = mi{N\uu + dd\N), (2) 
g s =m s {N\ss\N), (3) 

where mi = {m u + rrid)/2. 

Following the technique described in Refs. [3 [20], we fit 
octet baryon mass data recently published by the PACS- 
CS Collaboration [5T] using a chiral expansion: 

M B = M {0} + 5M B 1] + 5M ( g /2) + ... (4) 

Here, denotes the degenerate mass of the baryon 

octet in the SU(3) chiral limit, 5M B ^ gives the correc- 

(3/2) 

tion linear in the quark masses, and 5AI B represents 
quantum corrections corresponding to one-loop contri- 
butions from the pseudo-Goldstone bosons <j) — 7r, K, rj. 
Explicit expressions for the extrapolation formulae and 
for renormalisation of the loop integrals may be found in 
Ref. [D]. 

Following Ref. [5], we retain the octet-decuplet mass 
difference 5 in numerical evaluations to properly account 
for the branch structure near ~ S. The loop contribu- 
tion parameters are set to appropriate experimental and 
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phenomenological values; D + F = gA = 1.27, F = l-D, 
C* = -2D, f = 0.0871 GeV, and 5 = 0.292 GeV. Within 
the framework of FRR, we introduce a mass scale A, 
through a regulator u(k). A is related to the scale beyond 
which a formal expansion in powers of the Goldstone bo- 
son masses breaks down. In practice, A is chosen by fit- 
ting to the lattice data itself. For further discussions of 
the FRR regularization scheme, we refer to Refs. [22Tl26j . 
To provide an estimate of the model-dependent uncer- 
tainty in our result, we consider a variety of forms of the 
regulator u(k), namely monopole, dipole, and Gaussian, 
as well as a sharp cutoff. To further estimate systematic 
uncertainties, we allow /, the meson decay constant in 
the chiral limit, the baryon-baryon- meson coupling con- 
stants F and C, and 6 to vary by ±10% from the central 
values given above; see Ref. [27] for details. The effect of 
these variations are included in the final quoted errors. 

The PACS-CS results have been corrected for small, 
model-independent, finite volume effects before fitting. 
These finite volume corrections were evaluated by con- 
sidering the leading one-loop results of chiral EFT [SJ I251 - 
[30]. We note that the largest shift was -0.022 ± 
0.002 GeV for the nucleon at the lightest pion mass. 

The fit to the PACS-CS baryon octet data is shown in 
Figure [l] We find an optimal dipole regularization scale 
of A = 0.9 ± 0.1 GeV, in close agreement with the value 
deduced from an analysis of nucleon magnetic moment 
data [31] and, from the phenomenological point of view, 
remarkably close to the value preferred from comparison 
of the nucleon's axial and induced pseudoscalar form fac- 
tors [32]. The minimum xlof is °- 41 ( 6 -l/(20-5)) for the 
dipole, and varies between 0.40 and 0.42 for the other 
regulators. This value is somewhat lower than unity, as 
correlations between the lattice data cannot be accounted 
for without access to the original data. 

Clearly, the fit is very satisfactory over the entire range 
of quark masses explored in the simulations. Further- 
more, the masses of the octet baryons agree remarkably 
with experiment at the physical point. A comparison 
of the extrapolated baryon masses with the best experi- 
mental values is given in Table [TJ The first error quoted 
is statistical and includes the correlated uncertainty of 
all of the fit parameters including the regulator mass A, 
while the second is an estimate of model-dependence. 
This includes the full variation over dipole, monopole, 
sharp cutoff and Gaussian regulator forms, as well as ac- 
counting for the variation of the phenomenologically-set 
parameters F, C and S described earlier. 

As we fit baryon mass functions to lattice data over 
a range of pseudoscalar masses significantly larger than 
the physical values, it is prudent to check the consistency 
of our results as the analysis moves outside the power- 
counting regime (PCR), where higher order terms may 
become significant. By performing our fit to progressively 
fewer data points, that is, by dropping the heaviest mass 
points, we test the scheme dependence of our evaluation. 
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FIG. 1. Fit to the PACS-CS baryon octet data. Error bands 
shown are purely statistical, and incorporate correlated un- 
certainties between all fit parameters. Note that the data 
shown has been corrected for finite volume and the simula- 
tion strange quark mass, which was somewhat larger than the 
physical value. The green stars show experimental values. 

The results are consistent, and largely independent of 
the truncation of the data. This can be seen clearly in 
Figure [2j which shows the variation of the dimensionless 
baryon sigma terms as progressively fewer data points are 
used for the fit to the octet masses. The points shown 
correspond to an evaluation with a dipole regulator, and 
error bars are purely statistical. 
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Experimental 
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OBs 
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0.959(24) (9) 


0.939 


0.047(6) (5) 


0.022(6)(0) 


A 


1.129(15)(6) 


1.116 


0.026(3)(2) 


0.141 (8) (1) 


£ 


1.188(H)(6) 


1.193 


0.020(2)(2) 


0.172(8)(1) 




1.325(6) (2) 


1.318 


0.0089(7) (4) 


0.239(8)(1) 



TABLE I. Extracted masses and sigma terms for the physical 
baryons. The first uncertainty quoted is statistical, while the 
second results from the variation of various chiral parameters 
and the form of the UV regulator as described in the text. 
The experimental masses are shown for comparison. 



To further test our claim that the fitted mass functions 
accurately describe the variation of the baryon masses 
with quark mass, we compare our extrapolation with in- 
dependent lattice data along a very different trajectory in 
the mi — m s plane, as compared to the fit domain. Most 
lattice simulations, including that of the PACS-CS Col- 
laboration, hold the simulation strange quark mass fixed 
near the physical value, and progressively lower the light 
quark mass to approach the physical point. However, the 
UKQCD-QCDSF Collaboration has recently presented 
an alternative method of tuning the quark masses, in 
which the singlet mass (2m|- + m%) is held fixed [33] . 
This procedure constrains the simulation kaon mass to 
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FIG. 2. Dimensionless baryon sigma terms, evaluated using a 
dipole regulator, based on fits to the NPLQCD results at the 
lightest 5 (all), 4, and 3 pseudoscalar mass points. 
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FIG. 4. Locations of lattice QCD simulations by the PACS- 
CS Collaboration (blue circles), and UKQCD-QCDSF Collab- 
oration (red and green squares and diamonds) in the mi — m 3 
plane. The star denotes the physical point. Figure [3] shows 
the fit to the PACS-CS data only, evaluated at the UKQCD- 
QCDSF simulation quark masses. 



1.20 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 

2 / v 2 



FIG. 3. Prediction of UKQCD-QCDSF lattice data, based 
on our fit to the PACS-CS octet baryon mass simulation. 
Red (square) and green (diamond) points correspond to 24 s 
and 32 3 lattice volumes respectively. Error bands shown are 
purely statistical, and incorporate correlated uncertainties be- 
tween all fit parameters. 

always be less than the physical value. In comparison, 
the traditional trajectory in the — m^ plane neces- 
sarily keeps the kaon mass larger than the physical value. 

The close match between our fit to the PACS-CS points 
and the UKQCD-QCDSF lattice data, shown in Figure|3j 
is extremely encouraging. We emphasize that the lines in 
Figure[3]are not a fit to the data shown, but rather a pre- 
diction, resulting from the described fit to the PACS-CS 
octet data being evaluated along the UKQCD-QCDSF 
simulation trajectory. 

All lattice points shown in Figure [3] have been shifted, 
by the procedure described for the PACS-CS data, to ac- 
count for finite- volume effects. We chose to use the lattice 



spacing a = 0.078 fm deduced by the UKQCD-QCDSF 
Collaboration. For further details of the UKQCD- 
QCDSF data set, and the normalizations A/v, we 
refer to Ref. [33] . 

Figure [4] illustrates the significance of the prediction 
shown in Figure [3j While our fit was made only to the 
PACS-CS data, it successfully reproduces the UKQCD- 
QCDSF lattice results, at points in the mi — m s plane 
which arc substantially different from the PACS-CS sim- 
ulation trajectory. This very strongly supports our claim 
that the sigma terms, which correspond to derivatives in 
the plane shown in Figure [4j are accurately determined 
by our fit. 

To extract the sigma commutators from our baryon 
mass functions, we use the Feynman-Hellman rela- 
tion [31], 

dM B 

and, as above, replace quark masses by meson masses 
squared: mi — > and m s — » (m 2 K — rn^/2). For the 

case of the nucleon, we recall the alternative conventional 
notation to quantify the strangeness content, namely the 
kaon sigma term 

ctkn = \{mi +m s ){N\ll + ss\N). (6) 

As in Equation [T] overlines indicate corresponding (di- 
mensionless) quantities normalized by the baryon mass. 

A direct measure of the magnitude of the strange quark 
content of the nucleon relative to its light quark content: 

2(N\ss\N) = 
V (N\uu + dd\N) m s ^n 1 
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can be trivially evaluated given the strange and light 
quark sigma terms. At the physical point, we find 
cr„ N = 45 ± 6 MeV, <j K n = 300 ± 40 MeV and a s = 
21 ± 6 MeV, corresponding to a y-value of 0.04 ± 0.01 
for mi/m s = 0.039(6) [35] ■ This analysis also constrains 
0Vjv — fo to be 1.64±0.53 MeV. The quoted errors include 
all systematic and model-dependent uncertainties com- 
bined in quadrature. Results for the other octet baryons 
are made explicit in Table [Tj 

An advantage of the method used here is that we can 
easily evaluate sigma terms from our fit at any pion 
or kaon mass. The QCDSF Collaboration has recently 
presented very precise direct calculations of ovjv and 
a s from lattice QCD, at light quark masses somewhat 
larger than the physical values [T5]. At (m„mf{) values 
of (281,547) MeV, the Collaboration quotes er^ DSF = 
106(H)(3) MeV, and <jQ CDSF = 12±^ MeV. This 
compares very favorably to the results of our fits at 
these particular pseudoscalar masses, namely ovjv = 
131(H)(5) MeV and a s = 16(5)(1) MeV. Once again, 
the two uncertainties correspond to the described evalu- 
ation of systematic and model-dependent uncertainties. 

The conclusion of our analysis is clear. By develop- 
ing closed-form functions for baryon mass as a function 
of quark mass based on a fit to PACS-CS Collaboration 
lattice data, we were able to determine precise baryonic 
sigma terms by simple differentiation. This method al- 
lows us to achieve small statistical and model-dependent 
uncertainties. Moreover, we find excellent agreement 
with recent direct lattice calculations of these values at 
unphysical pseudoscalar masses. Decidedly our most sig- 
nificant result is a very precise value for the strangeness 
nucleon sigma term, namely a s = 21 ± 6 MeV at the 
physical point. 
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